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We construct a 2D quantum spin model that realizes an Ising paramagnet with gapless edge 
modes protected by Ising symmetry. This model provides an example of a "symmetry-protected 
topological phase." We describe a simple physical construction that distinguishes this system from 
a conventional paramagnet: we couple the system to a Z2 gauge field and then show that the 7r-flux 
excitations have different braiding statistics from that of a usual paramagnet. In addition, we show 
that these braiding statistics directly imply the existence of protected edge modes. Finally, we 
analyze a particular microscopic model for the edge and derive a field theoretic description of the 
low energy excitations. We believe that the braiding statistics approach outlined in this paper can 
be generalized to a large class of symmetry-protected topological phases. 



I. INTRODUCTION 

Recently it was realized that there are two distinct 
types of time reversal invariant band insulators: topo- 
logical insulators and conventional insulators .1^^ The two 
families of insulators are distinguished by the fact that 
topological insulators have protected gapless boundary 
modes while trivial insulators do not. Importantly, both 
time reversal and charge conservation symmetry play a 
crucial role in this physics. If either of these symme- 
tries is broken (either explicitly or spontaneously), the 
boundary modes can be gapped out and the sharp dis- 
tinction between topological insulators and conventional 
insulators disappears. 

This observation motivates a generalization of topo- 
logical insulators called "symmetry-protected topologi- 
cal (SPT) phases'^f^^Wi^. Consider a general gapped quan- 
tum system. The system may be built out of fermions 
or bosons/spins, and can live in any spatial dimension. 
Suppose that the Hamiltonian is invariant under some set 
of internal (on-site) symmetries. We will say that such 
a system belongs to a nontrivial SPT phase if it satisfies 
three properties. The first property is that none of the 
symmetries are broken spontaneously - the ground state 
is featureless from a symmetry-breaking point of view. 
The second property is that the ground state belongs to 
a distinct quantum phase from a "trivial state" with the 
same symmetry. That is, one cannot continuously deform 
the Hamiltonian in such a way as to connect the ground 
state with a "trivial state" without breaking one of the 
symmetries or closing the energy gap. Here, by a "trivial 
state", we mean a product state (in the boson/spin case) 
or an atomic insulator (in the fermion case). The final 
property of an SPT state is that the difference between 
it and the trivial state disappears if any of the symme- 
tries are broken. More precisely, the ground state can 
be connected with a trivial state without closing the en- 
ergy gap if one or more of the symmetries are broken 
during the process. We note that SPT phases typically 
exhibit robust gapless boundary modes analogous to that 



of topological insulators, though we will not include this 
property in the formal definition. 

Symmetry-protected topological phases have a long 
history in the one dimensional (ID) case. Most fa- 
mously, the Haldane phase of the 5 = 1 Heis enberg 
antiferromagnelji^ is known to belong to this clas^^'^^ESI 
More recently, a complete classification of ID SPT 
phases was obtaine d for both boson/spin systems-SAIll 
and fermion systems.'^E^l 

Much less is known about higher dimensional SPT 
phases. In the case of fermion systems, our understand- 
ing is largely limited to non-interacting models such as 
topological insulators or superconductors. For these sys- 
tems, an (almost) complete classification of SPT phases 
was obtained by Ref. [17] and UHl In some cases, it is 
known that this classification scheme is not affected by 
interactions (e^ the Z2 classification of topological in- 
sulators in twol^ and thre^^MU dimensions) . In general, 
however, this need not be the cas^^l ^nd consequently 
our understanding of interacting fermionic SPT phases 
in higher dimensions is incomplete. 

The boson case has received even less attention, and 
will be our focus here. In this case, a major advance 
was made by the recent paper, Ref. (Ill In that paper, 
the authors proposed a general classification scheme for 
bosonic SPT phases in general spatial dimension. Also, 
the authors constructed concrete microscopic models re- 
alizing each of these phases. This work established that 
the boson case is tractable even for interacting systems. 

Nevertheless, a number of questions remain open. One 
problem is that we have not identified any physical 
properties that distinguish different SPT phases in the 
hulk. The boundary physics is also poorly understood: 
while Ref. [10] showed that the 2D SPT states have 
symmetry-protected gapless boundary modes, the prob- 
lem for higher dimensions remains open. 

In this work, we address these (and other) questions in 
the context of a simple example. Specifically, we consider 
the case of 2D spin systems with a Z2 Ising-like symme- 
try. According to Refs. [10] and [TTl there is exactly one 
nontrivial SPT phase with this symmetry. This phase 
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can be thought of as a new kind of Ising paramagnet. 
Here, we construct an exactly soluble spin model that re- 
alizes this phase. We then derive three main results. Our 
first result is a simple argument that this model belongs 
to a distinct phase from a conventional Ising paramag- 
net. We derive this result by coupling the model to a Z2 
gauge field. After following this procedure, we find that 
the resulting gauged spin model supports quasiparticle 
excitations with different braiding statistics from that of 
a conventional (gauged) paramagnet. More specifically, 
we find that in a conventional paramagnet, the tt-Aux 
excitations have bosonic or fermionic statistics, while in 
the new paramagnet they have semionic statistics. It 
then follows immediately that the two paramagnets can- 
not be continuously connected without breaking the Zi 
symmetry or closing the energy gap. Closely related to 
this observation, we show that the two spin models are 
"dual" to two previously studied lattice models - each of 
which realizes a different type of Z2 gauge theory. This 
duality establishes a connection between SPT phases and 
previous worlP^l on the classification of topological gauge 
theories. 




FIG. 1. The Hamiltonians Hq,Hi (T]2i for the two spin 
models, (a) The Hamiltonian Hg is a sum of single spin terms, 
CTp. (b) The Hamiltonian Hi is a sum of seven spin terms 



{pii') 



where the product runs over the 



six triangles {pqq') containing p. 



II. TWO KINDS OF ISING PARAMAGNETS 

To begin, consider the following spin-1/2 model defined 
on the triangular lattice (Fig. [T^): 



(1) 



Our second result is a proof that the new paramagnet 
has gapless edge modes protected by Ising symmetry. In- 
terestingly, our argument reveals that the protected edge 
states are deeply connected to the braiding statistics of 
the TT-fluxes. This approach to proving edge state protec- 
tion is somewhat different from the original argument of 
Ref. [inland may be more amenable to higher dimensional 
generalizations. In the final part of the paper, we ana- 
lyze the protected edge modes at a more concrete level, 
focusing on a particular microscopic model of the edge. 
We derive a field theoretic description of the low energy 
modes, and analyze their stability to perturbations. 



Although we focus our discussion on a particular SPT 
phase, we believe that our basic approach is more general. 
That is, we expect that in a large class of SPT phases, 
braiding statistics can be used to uniquely characterize 
the bulk and to derive the existence of protected bound- 
ary modes. We discuss these potential generalizations in 
the conclusion. 



This paper is organized as follows. In section |TTj we 
describe spin models that realize both the conventional 
and the new kind of Ising paramagnet. In section [ITT] we 
show that the two spin models can be distinguished by 
the braiding statistics of the tt-Aux excitations. In section 
|IV| we show that the two spin models are dual to two 
previously studied lattice models. In section |Vj we show 
that the 7r-flux braiding statistics are directly connected 
to the existence of protected edge modes. Finally, in 
section [Vl] we analyze a particular microscopic model for 
the edge. 



This model describes a (conventional) Ising paramagnet. 
To see this, note that the system satisfies two prop- 
erties. First, the Hamiltonian is invariant under the 
Ising symmetry 5" = rip'^^p- Second, the ground state 
l^'o) = \(7p — 1) is gapped and unique - implying that 
the symmetry is not broken spontaneously. 

Surprisingly, there is another type of Ising paramagnet 
which is qualitatively different from Hq and represents a 
distinct quantum phase. A microscopic model for this 
new type of paramagnet was first constructed in Ref. 1101 
Here we describe another model which is more conve- 
nient for our purposes. The model we consider is a spin- 
1/2 system on the triangular lattice. The Hamiltonian is 
given by (Fig. 



Hi 



{pqq') 



(2) 



where the product runs over the six triangles {pqq') con- 
taining the site p. We note that this Hamiltonian is Her- 
mitian despite the factors of i. To see this, notice that the 
product includes a factor of i for each pair of neighboring 
spins (7, q' that have opposite values of a^. In particular, 
since the number of such pairs is necessarily even, the 
product always reduces to a factor of ±1. It is then clear 
that hI = Hi. (For readers who are c urio us as to how 
this model was constructed, 
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see section 

First we show that Hi describes a paramagnetic phase 
- that is the Ising symmetry is not spontaneously bro- 
ken. To establish this fact, we solve Hi explicitly. For 
simplicity, consider a geometry with periodic boundary 
conditions (i.e. a torus) . The key point is that 
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FIG. 2. A schematic plot of the ground states 'I'o and 'I'l 
for the two paramagnets Ho,Hi. (a) In terms of domain 
wall configurations, the ground state 'to is a equal weight 
superposition of all configurations, (b) The ground state "I'l 
is also a superposition of all domain wall configurations, but 
each configuration enters with a sign l)'^<i™ where Ndw is 
the total number of domain walls. 




(a) (b) 

FIG. 3. The Hamiltonians Ho, Hi jsjl for the two gauged 
spin models, (a) The Hamiltonian Ho is a sum of two terms. 
The first term is the gauge flux term i-ipqjJ,qrlJ-rp (thick tri- 
angle) where fip^ denotes the Z2 gauge field on the link 
(pq). The second term is the spin interaction cjpOp where 
Op — Y[(pqr) (1 + Mp<jM9rMrp)/2 and the product runs over the 
six triangles adjacent to p. (b) The Hamiltonian Hi includes 
the same gauge flux term /ipq/iq^/i^„ but has a more compli- 
cated seven spin interaction BpOp Kt. 



III. COUPLING THE SPIN MODELS TO A Z2 
GAUGE FIELD 



as can be verified by straightforward algebra. As a result, 
we can simultaneously diagonalize {Bp}. We will label 
the simultaneous eigenstates by |{6p}) where bp — ±1 
denotes the eigenvalues of Bp . It is not hard to show that 
there is a unique state for each choice of {bp}. In other 
words, the {bp} are a complete set of quantum numbers. 
We therefore have the full energy spectrum: each state 
\bp) is an energy eigenstate with energy 

E = -J2bp (4) 

p 

In particular, we can see that the ground state = 
\bp = 1) is unique and gapped. These two properties 
then imply that the Ising symmetry is not spontaneously 
broken. 

It is illuminating to compare the ground state wave 
functions of Hq, Hi. The ground state of Hq is the state 
where ct^ = 1 everywhere. Working in the basis, the 
wave function is given by 

*o({ap}) = 1 (5) 

for all spin configurations {up —t,i} (Fig- [2^)- As for 
Hi, we note that the ground state is the unique state 
with Bp = 1 everywhere. It is straightforward to check 
that the corresponding wave function is given by 

vI/i(K}) = (-l)^- (6) 

where Ndw is the total number of domain walls in the 
spin-configuration {ap =tii} (Fig. [2]d). We can see that 
the two ground states are nearly identical, differing only 
by some phase factors. Nevertheless, these two states 
belong to two different quantum phases, as we now show. 



In this section, we show that Hq, Hi belong to distinct 
quantum phases. Our strategy is as follows. Because 
Hq , Hi have a Z2 symmetry, we can couple them to a 
Z2 gauge field fi^^ = ±1 which lives on the links {pq) of 
the triangular lattice. We then show that the resulting 
gauged spin models have quasiparticle excitations with 
different braiding statistics. More specifically, we show 
that the two systems differ in the statistics of the tt-Aux 
excitations: while the 7r-fluxes have bosonic or fermionic 
statistics in the case of Hq, they have semionic statis- 
tics in the case of Hi. It then follows immediately that 
Hq , Hi cannot be continuously connected without break- 
ing the Z2 symmetry or closing the energy gap. 

Coupling Hq ,Hi to a Z2 gauge field requires several 
steps.l^ The first step is to apply the minimal coupling 
procedure, replacing nearest neighbor spin-spin interac- 
tions like (JqCTq, with a-qfJ-qq/O-q, . Next, we multiply each 
term in the resulting Hamiltonian (either ap or Bp) by 
the operator 

(pgr) 

where the product runs over the six triangles (pgr) ad- 
jacent to site p. The operator Op is a projector which 
projects onto states that have vanishing fiux through each 
of the adjoining triangles. We include this projection op- 
erator in order to ensure that our gauged Hamiltonian is 
Hermitian, and also to make the minimal coupling proce- 
dure unambiguous. (For more general models, we would 
replace Op with an operator that projects onto states 
that have vanishing flux through all the triangles in the 
vicinity of the spin-spin interactions) . The final step is to 
add a term of the form - J2{pqr) t^pqf^qrt^rp to the Hamil- 
tonian. This term ensures that the states with vanishing 
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^2 flux have the lowest energy. The resulting models are 
given by (Fig. [3|: 



P (pqr) 
i/l = - ^ BpOp - ^ nlqU'^qr^J.lp 



{pqr) 



where 



Bp = -o; n * — ^ 



(8) 



(9) 



Like all gauge theories, these models are defined on a 
Hilbert space consisting of gauge invariant states - that 
is, all states satisfying the constraint 



n 



PI 



(10) 



for all sites This constraint can be thought of as a 
Z2 analog of Gauss' law, V ■ E — Airp. 

Importantly, all the terms in _ffo, Hi commute with one 
another so these Hamiltonians can be solved exactly just 
like the ungauged spin models _ffo. Hi. In particular, it is 
easy to verify that both models have a finite energy gap. 

The next question is construct the quasiparticle excita- 
tions and show that they have different braiding statistics 
in the two systems. The quickest way to derive this fact is 
to note that Hq, Hi can be exact ly mapped onto the pre- 
viously studied "toric code"'2Sl2SI and "doubled semion"'^ 
models. These two models have been analyzed in detail 
and are known to support quasiparticle excitations with 
different statistics.'^SJ A description of these m odel s as well 
as the mapping to Hq , Hi is given in section 
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Alternatively, we can directly compute the quasiparti- 
cle statistics of Hq^ Hi and show that they are different. 
The first type of excitation is a "spin-flip", which we 
will denote by e. These excitations correspond to sites 
p where ct^ = — 1 for the case of _ffo, or Bp — —1 for 

the case of Hi. The second type of excitation is the 
"TT-flux" , m. These excitations correspond to triangular 
plaquettes (pgr) where ^^pq^J^qr^^rp = ^1- i&ct, there 
are two types of 7r-flux excitations, which differ by the 
addition of a spin-flip: to;, — nia ■ e. 

It is clear that in both systems, if we braid a spin- 
flip excitation e around either of the 7r-flux excitations 
TTia, mb, the resulting statistical Berry phase is tt (in some 
sense this is the definition of a 7r-flux excitation). It is 
also intuitively clear that the spin-flip excitation e is a 
boson in both models. All that remains is to understand 
the statistics of the 7r-fluxes. As we will now show, this 
is where the two models differ. 

To determine the 7r-flux statistics, we first identify op- 
erators that create these excitations. Like all quasipar- 
ticles with nontrivial braiding statistics, the vr-fluxcs can 
be created using an extended string-like operator. ■^^ If 
we apply these string-like operators to the ground state. 




FIG. 4. The string operator Vg (111 is defined for any path 
/3 on the dual honeycomb lattice and is given by a product of 
fj,pg over all links (pq) crossing (3 (thickened lines). Applying 
this operator to the ground state |^o) creates two vr-fluxes at 
the endpoints of /3 (shaded triangles). 



the result is a pair of 7r-flux excitations - one at each 
end of the string. In the case of Hq, the following string 
operator does the job: 



yfi - 



n 



(11) 



Here /? is a path in the dual honeycomb lattice joining 
the two triangular plaquettes, and the product runs over 
all links {pq) crossing /3 (Fig. |4]). We can verify that 
creates 7r-flux excitations at the two endpoints of (3 
by noting that Vg anticommutes with the flux ^p^^^^^^p 
through the two triangles at the ends of /3. At the same 
time, this operator commutes with all the other terms 
in Hq so it does not create any additional excitations.'^Sl 
Closely related to this fact, one can check that the state 
F^l^'o) does not depend on the choice of path 13, but only 
on the endpoints of /3 - a general feature of such string- 
like operators .'^SlsZ! \Ye will denote the 7r-flux excitation 
created by Vg by ma- A similar string operator creates 
the other type of 7r-flux, = i^a ■ e- 

In general, one of the most important aspects of string 
operators is the commutation relations satisfied by two 
intersecting strings. Let 7 be two paths on the dual 
honeycomb lattice that intersect one another. Using the 
definition (11), we can see that the two corresponding 



string operators commute with one another: 



P 7 7 



(12) 



This string algebra is important because we can use it 
to find the statistics of the quasiparticle nia (see Refs. 
[25H?7l for a general discussion of this connection) . One 
way to see this is to consider the special case where /3 is 
a closed path and 7 is an open path, as in Fig. [Sj In this 
case, the two operators Vp and have different physi- 
cal interpretations: while the operator can be thought 
of describing a physical process in which two 7r-fluxes are 
created and then moved to the endpoints of 7, the opera- 
tor Vp does not create any excitations at all. In fact, it is 
easy to check that exactly commutes with the Hamil- 

tonian Hq whenever (3 forms a closed loop. This suggests 
that should be thought of as describing a three step 
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FIG. 5. A schematic picture of the two states Vg V^l^fo), 
V^Vg l^'o). The first state (left) is obtained from a process 
in which two vr-fluxes are created at the endpoints of 7, and 
then two more fluxes are created, braided around the path /3 
and then annihilated. The second state (right) corresponds 
to executing these two steps in the opposite order. We expect 
these two states to differ by the Berry phase e^'" associated 
with braiding one 7r-flux excitation around another. The same 
is true for V/s , . 



process in which (1) two 7r-fluxes are created, (2) one of 
the TT-fluxes moves all the way around the closed path /3, 
and then (3) the two 7r-fluxes are annihilated. Using this 
interpretation, we can see that the state y^V^lvI/) is the 
end result of a process in which two 7r-fluxes are created 
at the endpoints of 7, and then afterwards another tt-Aux 
is braided around one of the endpoints and annihilated 
with its partner. In contrast, the state V^VjI^E') corre- 
sponds to executing these two steps in the opposite order. 
Comparing these two processes, we expect that they will 
differ by a phase factor which is exactly the statistical 
Berry phase associated with braiding one tt-Aux around 
another. In other words, the phase difference between 
these two states should be e^*^ where 9 is the exchange 
statistics for the particles: 



(13) 



In light of this relation, ( 12 ) implies that 9 — or ir. That 
is, nia is either a boson or a fermion. A similar analysis 
shows that the other tt-Aux excitation, rUb, is also either 
a boson or fermion. In fact, with a bit more work one 
can establish the more precise result that iria is boson 
and mf, is a fermion. The difference in statistics between 
Too, nib comes from the fact that rub = rua • e where e, rua 
have mutual statistics tt. However, we will not need this 
more detailed result here. (See Refs. 1^ and [57] for an 
analogous calculation for the closely related "toric code" 
model). 

We can repeat the same analysis for Hi. In this case, 
the following string operator creates a 7r-flux excitation: 



(14) 



yp= n n ^ 

{pq)i-P {pqq'),r 

■ n (-1^"'- n (i+m>,va^p,')/2 

{pqq'),i {pqq')ei3 



Here, the first product runs over all links (pq) crossing /3. 
The next two products run over all triangles (pqq') along 
the path such that g, q' are to the right of /3 or to the left 
of (3 respectively (Fig. [6]). The last product runs over all 




FIG. 6. (Color online) The string operator ( 14 1 is deflned 
for any path /3 on the dual honeycomb lattice. It acts on 
all triangles {pqq') along the path /3 (thickened lines). The 
action is different depending on whether q, q' are to the left 
of P (purple sites) or to the right of /3 (blue sites). Applying 
this operator to the ground state l^*!) creates two vr-fluxes at 
the endpoints of /3 (shaded triangles). 



triangles along /3. The operator Spqq' is defined by 



-(l-'7>p,f^|)(l + f^>p 



(15) 



As in the previous case, one can check anticommutes 
with the flux /Zp^/z^^/i^p through the two triangles at the 

ends of /3, but commutes with the Hamiltonian Hi every- 
where else. Hence, if we apply to the ground state, it 
creates 7r-fluxes at the two endpoints of /3. We will again 
denote this tt-Aux excitation by ma- 
in this case, one can check that the string operators 
satisfy a slightly diAerent algebra: for any two paths /?, 7 
intersecting one another, we have 



y p y^ y p 



(16) 



Therefore by the same reasoning as in (13), we conclude 
that the statistical angle 9 satisAes 29 = tt, so that 9 = 
±7r/2. In other words. Too is a semion. A similar analysis 
shows that the other tt-Aux excitation rub is also a semion. 
With a bit more worlP^, one can show that Too, vrib have 
opposite statistics - that is 9 = 7r/2 in one case and 
9 = — 7r/2 in the other - but again we do not need this 
more detailed result here. 

We have shown that the 7r-Auxes have diAerent statis- 
tics in the two gauged spin models: these excitations are 
bosons or fermions in the case of Hq, and are semions 
in the case of Hi. This result provides a simple physical 
distinction between the two systems. It also proves that 
the two spin models Hq, Hi cannot be continuously con- 
nected with one another without breaking the Z2 sym- 
metry or closing the energy gap. Indeed, if such a path 
existed, then we could construct a corresponding path 
connecting the gauged spin models Hq, Hi - a contradic- 
tion. We note, however, that the above argument does 
not rule out the possibility of connecting Hq,Hi if the 
Ising symmetry is broken during the process. Indeed, in 
appendix [A] we construct an explicit path H{s) of this 
kind. 
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IV. DUALITY BETWEEN SPIN MODELS AND 
STRING MODELS 

In this section we explain the relationship between the 
spin Hamiltonians Hq, and previously known models. 
Specifically, we show that -ffo ; Hi are related via a duality 
map to two previ ously studied lattice models - the "toric 
code" modepSl^lland the "doubled semion" model.l^The 
latter two models are sometimes called "string models" 
and are special cases of the general class of "string-net" 
models constructed in Ref. 26, This duality provides an- 
other point of view on the braiding statistics analysis in 
the previous section, and also suggests a natural classi- 
fication scheme for general 2D bosonic SPT phases with 
finite unitary symmetry groups. 

We begin by defining the duality map: we note that 
every spin configuration {ct^ = ±1} on the triangular 
lattice defines a corresponding domain wall configuration 
on the honeycomb lattice. Formally, this correspondence 
is given by r^^ = (TpCTg where / is the link separating sites 
p, q and = =pl corresponds to the presence or absence 
of a domain wall. We will refer to these domain walls 
as "strings." An important point is that the dual string 
degrees of freedom always form closed loops - that is, 
they satisfy the condition Qy = 1 where (Fig. [?]) 



(17) 



lev 



Using this correspondence, we can map our spin Hamil- 
tonians Ho,Hi (l][2) onto dual string Hamiltonians: 




(18) 



/(ilcgs of ] 



These Hamiltonians are defined on a Hilbert space con- 
sisting of closed string states (i.e. states satisfying — 
1 everywhere). 

The dual Hamiltonians H^,Hf are closely related to 
two models studied in Ref. [211 the "toric code model"'2Sl 
(or more accurately, a modified version of the toric code 
model) and the "doubled semion" model. To understand 
the precise relationship, recall that the latter two models 
are defined on a Hilbert space consisting of all string 
states on the honeycomb lattice - both open and closed. 
The two Hamiltonians are (Fig. [t]) 

V p \iep I 



V p \iep /eiogs of p / 





FIG. 7. The toric code and doubled semion models Ht.c, Hd.s 
(191. In both systems, the Hilbert space is equivalent to a 
spin- 1/2 model where the spins live on the links / of the hon- 
eycomb lattice, (a) The toric code Hamiltonian Ht.c is a sum 
of two terms. The first term, Qv (17 1 is a product of over 
the three links adjacent to the vertex v. The second term 
involves the interaction Yli^ '''i which acts on the six links 
adjacent to the plaquette p. (b) The doubled semion Hamil- 
tonian Hd.a includes the same vertex term Qi,, but contains 
a more complicated plaquette term Higp t'i Hieicgs of p f^'^i) 
where ]{x) = i'^-^^/^ 



Here Pp denotes the projector Pp = HuspCl + Qv)l'2.. 
This operator defines a projection onto states that satisfy 
the closed string constraint Q^, = 1 at all vertices of the 
plaquette p. 

Comparing (19) and (18 1, we see that Hg,Hf can be 
obtained by restricting Ht.c,Hd.s to the closed string 
(Qv = 1) subspace. In other words, the spin models 
Hq , Hi are dual to a restricted variant of the toric code 



and doubled semion models. 

In fact, this duality can be extended to one that maps 
the gauged spin models Hq, Hi onto the unrestricted toric 
code and doubled semion models ( 19 1. The extended du- 
ality is defined by setting Tj 



^p^ql^pq.n 



fip where 



(19) 



I is the link separating sites p, q. Substituting these ex- 
pressions into Ht.c, Hd.s, it is easy to check that the result 
is exactly Hq,Hi. We note that this duality maps local 
operators onto local operators and should therefore be 
thought of as an exact equivalence between two quan- 
tum systems. Thus the gauged spin models Hq,Hi are 
physically identical to the toric code and doubled semion 
models. 

The above dualities are variants of the well-known cor- 
respondence between the 2D Ising model and 2D Z2 
gauge theory.l^lISl! Xo see this, note that the closed string 
models Hq, Hi are simply Z2 gauge theory Hamiltoni- 
ans, phrased in the language of strings. The Hamiltonian 
Hq is the standarcP^l22l Z2 gauge theory Hamiltonian (in 
the zero coupling limit where there is no electric energy 
term J^i while Hf is another kincF^of Z2 gauge the- 
ory. From this point of view, the correspondence between 
Hq, Hi and Hq, Hf is a duality between two types of 2D 
Ising paramagnets, and two types of 2D Z2 gauge theory. 

We can understand the duality between Hq,Hi and 
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fft.c, s in a similar way. We note that the first two 
models can be thought of as two types of Ising param- 
agnets coupled to (standard) Z2 gauge theory, while the 
latter two models can be thought of as two types of Z2 
gauge theory coupled to a (conventional) Ising paramag- 
net. Hence the duality between Hq,Hi and Ht,c,H(i.s is 
a variant of the well-known self-duality of 2D Z2 gauge 
theory coupled to Ising matter. 

We expect that these dualities can be generalized from 
Z2 to any finite unitary symmetry group G: each SPT 
phase with symmetry group G is dual to a correspond- 
ing gauge theory with gauge group G. This correspon- 
dence immediately suggests a classification scheme for 2D 
bosonic SPT phases with finite unitary symmetry groups: 
it is known that the different types of 2D gauge theo- 
ries with group G (or equivalently, different string-net 
models corresponding to G) are in one-to-one correspon- 
dence with elements of H^{G,U{1)). (For a derivation 
of this result, see Ref. '23;, and also section 10. I.E. 3 of 
Ref. 1301 ■ Hence, the duality map suggests that different 
SPT phases associated with symmetry group G can also 
be classified by H^{G,U{1)). This classification scheme 
is identical to the proposal of Ref. HU 

Another application of these dualities is that they give 
a simple method for constructing exactly soluble mod- 
els for bosonic SPT phases with finite unitary symme- 
try group G. The first step is to construct the differ- 
ent "string-net" model^^H corresponding to the group 
G. These are models with string types given by the 
group elements g d G, and branching rules given by 
group multiplication: {31,32, 33} is an allowed branch- 
ing if 31.9233 = 1- In general, there will be a finite 
number of different models with these branching rules 
- each one corresponding to a different solution of the 
self-consistency equations of Ref. |26l (For readers fa- 
miliar with the string-net models in Ref. 1261 we note 
that one can equivalently string-net models by letting the 
string types be irreducible representations of G, but these 
models are not as convenient for the duality construction 
when G is non-abelian). We then take the dual of these 
models, and thereby construct exactly soluble models for 
bosonic SPT phases. The models Ho,Hi discussed here 
were constructed using this approach. In appendix [B] 
we show that an analogous duality in a space-time La- 
grangian description can be used to construct topological 
non-linear sigma models for SPT phases. 



V. PROTECTED EDGE MODES AND 
BRAIDING STATISTICS 

The most dramatic distinction between the two types 
of paramagnets, Hq,Hi is that Hi has protected gap- 
less edge modes, while Hq does not. In other words, 
if we define Hi in a geometry with a boundary, then 
the energy spectrum always contains gapless excitations. 
These gapless excitations are guaranteed to be present 
as long as the Ising symmetry is not broken (explicitly 



or spontaneously). In this section, we give a general ar- 
gument proving this fact. Our argument reveals that 
these edge modes are closely connected to the semionic 
braiding statistics of the 7r-flux excitations in the gauged 
spin model. Hi . We note that the existence of protected 
edge modes was previously established in Ref. ITOl using 
a different approach. 



A. The argument 

The statement we prove is as follows. We consider a 
disk geometry with a Hamiltonian of the form 



H = H 



bulk 



H 



H 



bulk 



= -J2Bp (20) 



where Bp is defined as in ([2| and the sum runs over all 
sites p lying strictly in the interior of the disk. We take 
the edge Hamiltonian TJedgo to be any Ising symmetric 
Hamiltonian with local interactions which acts on the 
spins on or near the boundary of the disk. In this setup, 
it is clear that the ground state \'^) oi H satisfies Bp = 1 
when p is far from the edge; in fact, in order to simplify 
the discussion, we will assume that Bp ^ 1 for all p lying 
strictly in the interior of the disk. (This last assumption 
is not essential to the argument and is included only for 
expository purposes). Given these assumptions, we will 
show that 1^) cannot be both gapped and unique. We 
interpret this as evidence for a protected gapless edge 
mode. 

First, we present a formal mathematical argument 
proving this fact. Then, in section VB we discuss the 



physical interpretation of this argument as well as its re- 
lation to the braiding statistics of the 7r-flux excitations. 

Our argument is based on the following construction: 
let /3 be a path on the dual (honeycomb) lattice that 
joins two points a, b on the edge. We define an associated 
unitary operator Wp by 



wp^u-p- n ' 

p.int {pqq'),r 



n (-1) 

{p'iq'),i 



(21) 



Here, the first product runs over all sites p in the interior 
of the the path /3, while the last two products run over 
all triangles {pqq') along the path such that q, q' are to 
the right of /3 or to the left of (3 respectively (Fig. |8|. 
The operator Spqqi is defined by 



1 



opqq' 



(l-aX)(l + <^X') 



(22) 



(The reader may notice that the definition of Wp is simi- 
lar to that of the string operator Vp ( 14 ) . We discuss the 



relationship between the two operators in section V B ) 



We note that the unitary operator Wp satisfies several 
interesting properties: 

1. Wp is symmetric under S: S^'^WpS — Wp. 




a c b d 



a c b d 



FIG. 8. (Color online) The operator Wp (|2T| is defined for 
any path on the dual honeycomb lattice that joins points a, b 
on the edge. In the interior of the path /3 (shaded region), W/s 
acts like the symmetry transformation S — Yip <^p- The oper- 
ator also acts on all triangles {pqq') along the path P (thick- 
ened lines). The action is different depending on whether q, q' 
are to the left of /3 (purple sites) or to the right of P (blue 
sites) . 



2. In the (far) interior of (3, Wp acts like the symme- 
try transformation S — , while in the (far) 
exterior, Wp acts like the identity operator. 

3. Wp transforms local operators into local operators. 
That is, Wg^OWp is local if and only if O is local. 



4. The state W^^l^*) is independent of the choice of 
path p. 

5. Let O be a local operator which acts on spins within 
some convex region R not containing either of the 
endpoints of /3. Then O has the same expectation 
value in the two states |*) and |*') =Wp\'^). 

Properties [T] [2] are clear, while property [3] follows from 
the fact that Wp can be decomposed into a product of 
two sets of commuting local unitary operators. To derive 
property [4] we make use of a simple operator identity: for 
any two paths /3, /?' with the same endpoints. 



(23) 



where the product runs over sites p in between the two 
paths /3, /?'. Applying this identity, we deduce 

Wp\^)^Wp: Yl Bp\^)^Wp,\^) (24) 

pei3'-i3 

where the second equality follows from the fact that 
Bp 1 5') = for all p lying strictly in the interior of 
the disk. 

As for property [5] there are three cases to consider: 
the region of support R may be contained entirely in the 
exterior of /?, it may be contained entirely in the inte- 
rior, or it may overlap the path f3 itself. In the first 
case, Wp commutes with O, immediately implying the 
desired equality (*'|0|*') = {^\0\^). In the second 
case, Wp^OWp = S~^OS. Then, since l^*) is invariant 



y 



y 



FIG. 9. The string operators obey the anticommuta- 
tion relation WpW-y = —W-^/Wp for any two paths /3, 7 that 
intersect one another as above. 



under S, we again have {^'\0\^') = {^\0\^). The only 
case where the expectation value of O could be different 
in the two states is if R overlaps the path /3. However, 
by property [2] we can deform the position of /? with- 
out affecting W/jjvl/). In particular, we can deform f3 so 
that it avoids R. Therefore, the expectation values must 
coincide in this case as well. 

In addition to properties [l][5j it is easy to check that 
the Wp operators satisfy a nontrivial algebra: for any 
two paths /3 and 7 that intersect one another (see Fig. 
[9]) we have 



(25) 



We now show that this operator algebra, together 
with properties l][5 implies that cannot be a unique 
gapped ground state. The argument is a proof by con- 
tradiction. Suppose that j^*) is a unique gapped ground 
state. It then followJ^U that \^) has short-range corre- 
lations. That is, for any well separated local operators 
Oi,02, we have (^lOiOal*) = («'|0i|*)(«'|02|«') up to 
corrections which are exponentially small in the distance 
between Oi,02- Also |*') = must have short 

range correlations since Wp transforms local operators 
into local operators. 

Next, we use the fact that the two states j^*), W^l^*) 
share the same local expectation values away from the 
endpoints a, b (property [s]). Given this fact, and also the 
fact that both states have short range correlations, it fol- 
lows that we can find local operators Ua, C4 acting near 
a, h (or more accurately, exponentially localized opera- 
tors) such that 



(26) 



To see this, consider the analogous question for the con- 
ventional paramagnet |^'o) = |cr^ = 1): suppose that 
some state |5'q) with short-range correlations has the 
same local expectation values as jvPo) except near two 
points a, h. In this case, the state must have = 1 
far from a, 6, so it is clear that we can find local opera- 
tors Ua,Ub acting near a,b such that UaUbl'^'o) = l^'o)- 
Having established this property for |5'o), it follows that 
the same property must also hold for j^f) since j^*), |^'o) 
can be continuously connected without closing the en- 
ergy gap and are therefore equivalent up to a local uni- 
tary transformation.!^ (To see that the two states can 
be connected with one another, we note that appendix 
[a] gives an explicit path connecting j^*) and |4'o) in the 
bulk; as for the edge region, we make use of the general 
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any ID unique gapped spin state can be 
continuously connected with a ID product state). 

A key question is to understand how Ua, Ub transform 
under S. In appendix [C] we show that Ua, Ub can always 
be chosen so that they are either both even or both odd 
under S. Furthermore, this even or odd parity must be 
the same for all pairs of endpoints a, b. In other words, 
either all the Ux operators are even under S, or all of 
them are odd under S. 

We are now ready to complete the argument. In ad- 
dition to /3, we consider a second path 7 that connects 
two other points c, d on the edge. We choose these paths 
so that they intersect each other, and so that their end- 
points are well separated (see Fig. [9|. By the same 
reasoning as above, we have UcU^Wyl'^/) — j^*) for some 
local operators Uc, Ud acting near c, d. Now, define 



Wp = UaUbWp, W^^UJJdW^ 



By construction, 
Hence, 



I*) and W-,m 



(27) 



I*) 



(28) 



But at the same time, W jj, W ^ anti-commute, as we now 
show. To see this, recall that looks like S in the 
interior of /3 and the identity map in the exterior of /3 so 
that 

WpUMdW-^ = {SU,S-^)Ud = ±UcUd (29) 

where the sign is determined by the parity of Uc under 
S. Similarly, we have 



w^uMhW: 



UJSUbS-^) = ±UaUb (30) 



where the sign is determined by the parity of Ub under 
S. Importantly, these two signs are the same since the 
Ux operators all share the same parity. Hence, either the 
two pairs {Wfs^UcUd} and {W^,UaUb}, both commute or 
both anti-commute. In either case, the anti-commutation 



relation (25) implies that Wp,Wj anti-commute: 



(31) 



Comparing (281, (31 1, we arrive at a contradiction. 



Hence our assumption must be false and H cannot have 
a unique, gapped ground state. 



operator ( 14 ) constructed in section III 



Indeed the 



two operators are identical except for the fact that V^^ 
is written in terms of the formalism of the gauged spin 
model, while Wp is written in terms of the original un- 
gauged spin model. This similarity, together with the 
fact that creates 7r-fluxes at the two endpoints of /3, 
immediately suggests the above interpretation. 

Many of the properties of W13 can be understood using 
this interpretation. For example, property |4] is natural 
since we do not expect that the final state M^^j^') will 
depend on the path /? along which the 7r-fluxes are moved. 
Similarly, property |5] is sensible since the only effect of 
Wp is to create two 7r-flux excitations at the endpoints 
a, b. 

We can also understand the meaning of Wp — 
UaUbWp: this operator describes a process in which 
two TT-fluxes are created, moved to the boundary, and 
then subsequently annihilated. The annihilation at the 
boundary may require either an even or odd number of 
spin flips, corresponding to the two cases where Ua,Ub 
are even or odd under the symmetry S. Depending on 
whether the number of spin flips is even or odd, the op- 
erator Wp describes the annihilation of one or the other 
of the two types of 7r-fluxes. 

Most importantly, we can now see that the anti- 
commutation relation (31) is a direct reflection of the 



semionic braiding statistics of the 7r-flux excitations, just 
like (16). Given this interpretation it is easy to see how 
to generalize the proof to other systems. The only re- 
quirement for the above argument is that neither of the 
TT-fluxes are bosons or fermions: as long as this is the 
case, we can derive a nontrivial commutation relation 



2ie 



(32) 



where 6 is the exchange statistics of the 7r-flux. We can 
then deduce the existence of protected edge modes just 
as above. On the other hand, if one of the 7r-fluxes is a 
boson or a fermion - as in the conventional paramagnet 
Hq - the argument breaks down completely: the relation 
( 32 ) is no longer in contradiction with ( 28 ) . From this 



point of view, the key reason that Hi has a protected 
edge mode and Hq does not, is the difference in their 
TT-flux braiding statistics. 



VI. MICROSCOPIC EDGE ANALYSIS 



B. Physical interpretation 

In this section, we explain the physical meaning of the 
above argument and its relation to the braiding statistics 
of the TT-flux excitations. We begin with the heart of the 
argument - the unitary operator Wp. This operator has 
a simple physical interpretation: it describes a process 
in which two 7r-fluxes are created in the bulk and then 
moved along the path /? to points a, b at the boundary. 
To see this, notice that Wp is very similar to the string 



In this section, we investigate the protected edge 
modes of Hi at a more concrete level. We analyze a par- 
ticular example of a gap less edge for Hi, derive a field 
theoretic description of the low energy modes, and in- 
vestigate the effect of perturbations. As in section |Vj 
we consider a disk geometry, with a Hamiltonian of the 
form H — -ffbuik + Hedge- The bulk Hamiltonian i?buik is 
defined by i/buik = — J2p where the sum runs over all 
sites that are strictly in the interior of the disk. The edge 
Hamiltonian -ffcdgc can be any Ising symmetric Hamilto- 
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a„ a, a 





(a) 



FIG. 10. A schematic picture of the zero energy edge states in 
the case Hedge = 0. (a) The 2^ zero energy edge states can be 
parameterized by boundary spin configurations {ai, ...,ajv}, 
where a„ =t,4- (b) For each choice of {a}, the corresponding 
wave function 5'{a} is defined by '^{a}{{ciint}) = (—1)^''"' 
where Ndw is the total number of domain waUs in the system. 
We use a convention where we close up all the domain walls 
by assuming there is a "ghost" spin in the exterior of the disk, 
pointing in the f direction. 



nian with local interactions which acts on the spins on or 
near the boundary. 



in the exterior of the disk, pointing in the f direction 
(Fig. IOd). We will denote these states by |ai, aAr). 



As is apparent from this parameterization, we can think 
of these degenerate ground states as zero energy edge 
states. 

It is useful to define operators that act on 

\ai, un) just like the usual Pauli spin operators. We 
note that the (f*j operators should not be confused with 
the physical boundary spin operators aj, which act on the 
full Hilbert space of the spin system. In the case, the 
two types of operators are closely related - for example, 
(T^ = Poa^Po where Pq is the projection operator onto 
the 2^ dimensional edge state subspace. However, this 
simple relation does not hold for the or operators, 
or for more complicated products of spin operators. 

An important question is to understand how the sym- 
metry S acts on the edge states. Using the definition 
(33), one finds that the Ising symmetry S acts as 



N 



5|ai,...,ajv) = ± n '^L^J/^i'-'/^A') 



(34) 



where the sign depends on the configuration of q;„ i 
as follows: the sign is — if the total number of domain 
walls between the a's is divisible by 4 and + otherwise. 
In other words, the action of the Ising symmetry on the 
above basis states is described by the operator 



A. Zero energy edge states 

We begin with the case where H^dge — - that is, the 
edge Hamiltonian vanishes. In this case, we can compute 
the energy spectrum in the same way as we did for the 
periodic (torus) geometry. First, we simultaneously diag- 
onalize the Bp operators for all sites p that are strictly in 
the interior of the disk. Next, we note that each of these 
simultaneous eigenstates is an energy eigenstate with en- 
ergy E = — hp where &p = ±1 is the eigenvalue under 
Bp. The final step is to determine the degeneracy of these 
simultaneous eigenspaces. A natural guess, based on di- 
mension counting, is that each simultaneous eigenspace 
{hp = ±1} has a degeneracy of 2^, where N is the num- 
ber of spins along the boundary of the disk. In particular, 
we expect that there are 2^ degenerate ground states. 

We can verify this counting by constructing ex- 
plicit wave functions for these degenerate ground states. 
Specifically, we define a wave function ^'{q} for each 
boundary spin configuration {ai, a^}, where =t, i 
(Fig. 10 1). This wave function is a function of the spins 



dint i lying strictly in the interior of the disk, and is 
given by 



(33) 



where Ndw is the the total number of domain walls in the 
system. Here, we define Ndw using a particular conven- 
tion where we close up all the domain walls that end at 
the boundary by assuming that there is a "ghost" spin 



N 

n 



N 



■ exp 



n'-'n+ll 



(35) 



In order to gain some intuition about S, we note that 
the operators (7^,(t^,(7^ transform under the symmetry 
according to 



S = 
S-^alS = u 



' 71-1" n+l 



Ji" n+1 



(36) 



B. An example of an edge Hamiltonian 

We now imagine adding a nonvanishing edge Hamil- 
tonian Hedge- If ^^cdgc is small, then we can analyze its 
effect using degenerate perturbation theory. The first 
order splitting of the 2^ degenerate ground states can 
be obtained by diagonalizing Po^odgo-Po where Pq is the 
projection onto the zero energy edge state subspace. In 
general, PoHcdgcPa can be expressed as a function of the 
CT^ operators. We can therefore find the edge state spec- 
trum by so lvin g a ID spin chain with an unusual Ising 
symmetry ( 35 1 



Here we will focus on a particular choice of Hedge which 
can be solved exactly. We emphasize that this choice is 
not unique, and that other edge Hamiltonians may give 
rise to different edge spectra. Nevertheless, we believe 
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We conclude that: 



FIG. 11. The operator is defined just like Bp ([2]), except 
with an additional "ghost" spin in the exterior of the disk 
pointing in the f direction (dotted arrow). More explicitly, 
B^ acts on the three triangles ingq ') c ontaining the boundary 
spin n with an action given by (381. The operator B^^ is 
defined similarly. 



that the particular Hedge we consider is a useful illustra- 
tive example. We will derive a low energy edge theory for 
this example, and investigate the effect of perturbations. 
Specifically, we consider an edge Hamiltonian 



edge 



= - J 



N 

E 

n=l 



(37) 



where is defined just like Bp, except with an addi- 
tional "ghost" spin in the exterior of the disk, pointing 
in the ^ direction. More explicitly. 



n 

{nqq') 



(38) 



where the product runs over the three triangles {nqq') 
containing the boundary spin n, and where n -I- 1, n — 1 



denote the two neighboring boundary spins (Fig. 111. 
The operator B^ is defined the same way, except that we 
take the "ghost" spin to point in the J, direction. That 



B^ = -cr"" • ' 



n 

{nqq') 



(39) 



The above edge Hamiltonian has several nice prop- 
erties. First, the edge Hamiltonian is Ising symmetric. 
Indeed, this follows from the fact that S~^B^S — B^, 
S~^Bj^S = B^. Another property of ffedgo is that 
[iJedge, ^^buik] = 0. This property follows from the fact 
that the Bp operators commute with Bj^ , B^ (which in 
turn follows from the fact that the Bp operators commute 
with each other) . One consequence of this commutation 
relation is that the low energy edge spectrum obtained by 
diagonalizing Po^cdgc-Po is exact, rather than just being 
correct to first order in perturbation theory. 

We now compute Po^^odgc^o- Using the definition of 
the basis states ( 33 ) , it is easy to check that 



from which it follows that 



n-l" n" n+l 



(41) 



N 



PM 



O^cdgcJ^O 



n=l 



n" n+lJ 



(42) 



A nice feature of this Hamiltonian is that it has a U{1) 
symmetry: it conserves J2n=i ^nPn+i ~ the total number 
of domain walls between the boundary spins. In order 
to make this U{1) symmetry manifest and to simplify 
the analysis, it is useful to rewrite the Hamiltonian in 
terms of the dual domain wall variables. Naively, we can 
accomplish this by defining 



cr„tT 



n'-' n+l 



(43) 



and re-expressing everything in terms of the r's. How- 
ever, the above duality transformation doesn't quite work 
for a system with periodic boundary conditions, since 
the variables obey the global constraint n^^=i = 1; 
and therefore only describe — 1 independent degrees of 
freedom. (Equivalently, there is no way to express in 
terms of the variables) . 

In order to incorporate the missing degree of freedom 
and make the dual description complete, we introduce an 
additional Z2 gauge field n that lives on the links 
((n — l)n) connecting neighboring boundary sites (n — 
1), n. We then define the duality transformation between 
a and r, /i by the relation 



Mr, 



together with the gauge invariance constraint 



X X Z 1 

r-n—l.nh^n.n+l'n ^ 



(44) 



(45) 



It is easy to check that there is a one-to-one correspon- 
dence between configurations of (f ^ = ±1 and configura- 
tions of /ifj = ±l,r,j = ±1 obeying the constraint (45 1. 
Similarly, there is a one-to-one correspondence between 
physical operators written in terms of the (f's and gauge 
invariant combinations of /i, r (i.e. operators that com- 
mute with the left hand side of (45)). In particular, the 
operators ct'^,ct^,ct^ are given by 

'^n ~ ''"n-l'''nMn-l,n 

n n — 1 nPn— 1,71 

K = Mn-l,n (46) 

while the symmetry transformation S is given by 



N 



ri=l 



(47) 



71=1 



Using ( 46 1 to re-express the Hamiltonian ( 42 ) in terms 



(40) of the domain wall variables r, we find 



N 



Po^cdgcPo = -2J ^(r+„iT-Ai:_i_„ + h.c.) (48) 
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This Hamiltonian is the usual spin- 1/2 XX chain, cou- 
pled to a Zi gauge field /i^_i^„. The only effect of the 
Z-2 gauge field is to double the size of the Hilbert space 
so that is includes sectors with both periodic and anti- 
periodic boundary conditions for the r variables. The 
two types of boundary conditions correspond to the two 
possibilities n^=i A^fi-i.n = ±1- 



C. Edge theory 

Given previous work on the spin-1/2 XX chain, it is 
now straightforward to construct a field theory descrip- 
tion of the low energy edge modes. (We could also derive 
the exact edge spectrum, but this is less useful to us, 
as we ultimately want to analyze the effect of perturba- 
tions) . To be specific, the low energy excitations of the 
spin-1/2 XX chain (48) are knowrpS to be described by 
the non-chiral Luttinger liquid 



Combining these two results, we see that our expres- 



47r 



Stt 



(49) 



with Luttinger parameter K = 1, and velocity v — 4 J a 
where a is the lattice spacing. Here, we are using a nor- 
malization convention in which expressions of the form 
^ike+ii(p .^yj^jj integer k, I correspond to local spin oper- 
ators (i.e. gauge invariant combinations of T,fi). For 
example. 



'ri — 1 'n fn— l,n 



2a 



(50) 



In the above normalization convention, the boundary 
condition for 9 is that 9{L) = 9{0) (mod 2tt). This con- 
dition automatically incorporates both the periodic and 
anti-periodic sectors for r: the two sectors correspond to 
the two cases 9{L) ^ 9{0)+4mTT and 9{L) = 6'(0)-|-(4m + 
2)tt, as one can see using the heuristic t+ ~ e'*/^. The 
boundary condition for cj) is also 4>{L) = 0(0 ) (mod 27r). 

The last component of the edge theory ( [49| is to under- 
stand how the Ising symmetry S acts in the new variables 
9,(j>. To this end, we note that (501 implies that 



N 



\ n=l , 

Similarly, we have: 

N 



exp 1 / dr,(l)dx ] (51) 



n—1 



(52) 



This equality follows from the observation that the 
periodic/anti-periodic sectors 11^=1 Mn-i,n — ±1 cor- 
respond to the two boundary conditions 9{L) — 9{0) — 
4m7r, (4m -I- 2)tt respectively. 



sion (47) for 5* becomes 



S = exp ( ^ / dx9dx - ^ / dx0ax 



(53) 



(up to a phase factor). Using the commutation relation 
[9{x),dx4'{y)] = 2TTiS{x — y), we deduce that 



(54) 



The transformation law ( 54 ) together with the action 



( 49 ) gives a complete description of the low energy edge 



physics. 



D. Stability and instability of the edge modes 

In this section, we investigate the effect of perturba- 
tions on the gapless edge (49). We find two results. Our 



first result is that the edge is unstable: the edge modes 
can be gapped out by arbitrarily small Ising symmet- 
ric perturbations. Our second result is that the edge 
is protected: we find that in all cases where perturba- 
tions gap out the edge, the Ising symmetry is broken 
spontaneously In other words, we do not find any per- 
turbations which gap out the edge without breaking the 
Ising symmetry, explicitly or spontaneously. This result 
is consistent with the general edge protection argument 
presented in section |Vj 

We focus on a particular class of perturbations of the 
form 



U{li,l2) = U{x) cos{h9 + l2(t> - a{x)) 



(55) 



where /i,^2 are integers. Using (54), we can see that 
the perturbation U{li,l2) is even or odd under the Ising 
symmetry depending on whether li -\- 12 is even or odd, 
respectively. 

The above perturbations are all "local" in the sense 
that they can be generated by adding appropriate short 
range spin interactions at the edge. For example, the case 
t/(0, 1) can be generated by adding to the edge Hamil- 
tonian (37) a term of the form Ua^. Similarly, ?7(1,0) 
can be generated by the term U{Bl^ — B^). Higher val- 
ues of li , I2 can be generated by more complicated spin 
interactions. 

Of particular interest are perturbations of the form 
U{1, 0) and [/(O, I). We know from the standard analysis 
of the sine-Gordon model that these terms can drive the 
edge to a gapped state by freezing the value of 9 or (j). 
This gapping can occur even for infinitesimal ?7, if U{1^ 0) 
or C/(0, /) are relevant in the renormalization group sense. 

To determine whether any of these operators are rel- 
evant, we note that their scaling dimensions are given 
by: 



/i(/,0) 



h{0,l) 



(56) 



13 



Clearly, the smaller the value of I, the more relevant the 
perturbation. On the other hand, Ising symmetry (54) 



requires even I. Thus, the most relevant Ising-symmetric 
operators are C/(2, 0) and U{0,2). Setting K = 1, we 
see that the term U{2, 0) has a scaling dimension greater 
than 2 and is therefore irrelevant, but U{0, 2) is relevant. 
Hence, U{0, 2) describes an Ising-symmetric instability of 
the edge. 

Microscopically, the term U{0, 2) can be generated by 
adding a staggered spin interaction 



N 

uTx- 



(57) 



to Hedge (37). In this case, the resulting gapping of the 
edge modes can be analyzed exactly. The analysis is 
similar to the derivation above: first, one maps the per- 
turbed Hamiltonian onto an XX chain with a staggered 
coupling constant J„ = J+(— 1)"C/. Then, one solves the 
resulting system using a Jordan- Wigner transformation. 
One can check that the effect of the perturbation is to 
induce backscattering for the non-interacting fermions, 
and hence open up a gap of order U. 

Before proceeding further, we make two observations 
about this edge instability. The first observation is that 
the instability described by U{0, 2) requires the break- 
ing of discrete translational symmetry. Indeed, as dis- 
cussed in the previous paragraph, U{0, 2) corresponds to 
backscattering between the left and right moving Jordan- 
Wigner fermions. This backscattering process doesn't 
conserve the lattice momentum and therefore requires 
the breaking of discrete translational symmetry and the 
doubling of the unit cell. 

The second observation is that the edge instability per- 
sists for any value of the Luttinger parameter, K. To 
see this, note that ([56| implies that h{2, 0) • h{0, 2) = 4. 
There are three cases to consider: either (a) h(2,0) < 
2 < /i(0,2), (b) h{0,2) < 2 < h{2,0), or (c) h{2,0) ^ 
h{0,2) = 2. In the first two cases, either U{2,0) or 
U{0, 2) is relevant, implying that the edge is unstable. 
In the third case, both operators are marginal, but the 
edge is still unstable since small perturbations can affect 
K and therefore make either [7(2, 0) or C/(0, 2) relevant. 
This analysis implies that the above edge theory has an 
Ising symmetric instability for any value of K. (This 
instability is closely related to the fact that there is no 
stable algebraic long range ordered phase in the 2D Z4 
clock model"). 

Although the perturbation U{0, 2) can open up a gap 
at the edge, it also spontaneously breaks the Ising sym- 
metry. To see this, note that U{0, 2) drives the edge into 
a state where (j) is frozen at some fixed value. In such a 
state, the operator cos{(t> — a{x)) acquires a nonvanishing 
expectation value. But this operator is odd under S* (54) 



implying that the resulting state spontaneously breaks 
the Ising symmetry. This result is consistent with the 
general argument in section |V} the edge modes can never 
be gapped out without breaking the Ising symmetry, ei- 
ther explicitly or spontaneously. 



We have seen that the above edge is unstable in the 
sense that small perturbations can gap out the edge 
while simultaneously breaking the Ising symmetry. We 
do not know whether a different choice of edge Hamilto- 
nian Hedge can give rise to a stable edge. Nevertheless, 
whether or not a stable Z2 edge is possible, we believe 
that the generalizations of the Ising paramagnet Hi 
support stable gapless edge modes for n > 2. Our expec- 
tation is based on the following conjecture: we believe 
that the Z„ generalizations of Hi support edge modes 
described by ( 49 ) with a Z„ symmetry given by 



S-'^eS ^ e + 2n/n; S"^ V-^ = + 27rfc/n (58) 

with A: = 1, n — 1. In this scenario, the most relevant 
Zn symmetric perturbations are U{n,0),U{0,n). Exam- 
ining (56), we can see that both of these operators are 
irrelevant over the finite range 8/n^ < K < n? /2. Hence, 
if K lies in this range, then the edge is stable to small 
perturbations. 



VII. CONCLUSION 

In this paper we investigated a 2D bosonic SPT phase 
with a Z2 Ising-like symmetry. This SPT phase can 
be thought of as a new kind of Ising paramagnet. We 
showed that this phase can be distinguished from a con- 
ventional paramagnet by coupling the system to a Z2 
gauge field and then analyzing the braiding statistics 
of the TT-flux excitations. We found that while the tt- 
fiuxes have bosonic or fermionic statistics in a conven- 
tional paramagnet, they have semionic statistics in the 
new kind of paramagnet. This result immediately im- 
plies that the two types of paramagnets belong to distinct 
phases. We also showed that these semionic braiding 
statistics directly imply the existence of protected edge 
modes. To complete the picture, we analyzed a particu- 
lar microscopic edge model for this phase, derived a field 
theoretic description of the edge modes, and investigated 
their stability to perturbations. 

While this paper has focused on a particular example, 
we believe that our basic approach can be applied more 
broadly. The simplest extension would be to consider 
2D bosonic SPT phases with arbitrary unitary symmetry 
groups G. Following the same approach as above, we can 
couple each such SPT phase to a gauge field with gauge 
group G. We can then analyze the quasiparticle excita- 
tions in this system and find their braiding statistics. By 
analogy with the Z2 case, we expect that these braid- 
ing statistics can be used to uniquely characterize each 
SPT phase and to derive the existence of protected edge 
modes. The same approach could potentially be used for 
2D fermionic SPT phases with unitary symmetries. 

One can also imagine a generalization to higher dimen- 
sional bosonic/fermionic SPT phases with unitary sym- 
metries. Again, we envision coupling each phase to the 
appropriate gauge field and then analyzing the braiding 
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statistics in the resulting system. In the 3D case for ex- 
ample, we expect that gauged SPT phases will contain 
both particle-like and loop-like excitations. The analog 
of braiding statistics is then the Berry phase associated 
with braiding a particle-like excitation around a loop-like 
excitation. We find it plausible that these braiding statis- 
tics could be used to distinguish different SPT phases and 
to derive the existence of protected boundary modes just 
as in 2D. Similarly, the duality between bosonic 2D SPT 



phases and 2D gauge theories (section IV) may also ex- 
tend to higher dimensions. 

On the other hand, it is not clear how to apply these 
ideas to SPT phases with anti-unitary symmetries such 
as time reversal symmetry. The problem is that we do 
not know how to define a gauge field for an anti-unitary 
symmetry. This question, as well as the potential gener- 
alizations discussed above, is an interesting direction for 
future work. 
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Appendix A: Adiabatic equivalence of Ho, Hi in the 
absence of the symmetry 

While Hq, Hi cannot be continuously connected when 
the Ising symmetry is preserved, these two models can 
be connected when the symmetry is broken. Indeed, 
consider the one-parameter family of unitary transfor- 
mations 



Ue 



(pqr) 



(Al) 



where the product runs over all triangles (pgr). Us- 
ing these unitary transformations, we can define a one- 
parameter family of Hamiltonians H{6) = Ug^HoUo- 
These Hamiltonians have local spin-spin interactions and 
have a finite energy gap for any value of 9. Moreover, 
H{0) = Hq, and one can check that H{ti/2A) = Hi, us- 
ing the identity 



7r/24 



(A2) 



Hence, this construction gives an explicit path that con- 
nects Ho,Hi. We note that H{9) breaks the Ising sym- 
metry for intermediate values of 0, as required by the 



Appendix B: Topological non-linear sigma models 
for the two paramagnets 

In this section, we explore the duality between the 
spin models Hq,Hi and the string models Ht,c.-,Hd.s in 
a space-time Lagrangian formulation. Using this space- 
time duality, we construct topological non-linear sigma 
models describing each of the two paramagnet phases, 
and thereby make a connection to the analysis of Rcf. 

m 

We begin with the Lagrangian description of the toric 
code and doubled semion models, Ht.c, Hd.s- These mod- 
els - like all string- net model^^H- have a Euclidean space- 
time description in terms of Turaev-VircP^ invariants. 
In general, these invariants define a space-time partition 
function for any 3D manifold M and any triangulation 
of M into tetrahedra. For the above two models the 
Turaev-Viro invariants are of the form 



Z = 



1 

2^ 



ijk. 



n 



'-'kin ' 



(Bl) 



, link 



tetrahedron 



where the degrees of freedom i, j, k... live on the links of 
the tetrahedra and run over the finite set {0,1}. The 
variable denotes the number of vertices in the trian- 
gulation. 

To specify the two partition functions, we need to de- 
fine di, G]^l^ - the weights associated with the links and 
tetrahedra in the triangulation. For toric code model, we 
have do = c?i = 1 and^- 

^000 _ -1 
"-^000 ^ -'-I 

x-<000 _ ^110 _ ^011 _ /-<101 _ 1 
•-^111 — ^001 ^ "-^lOO ^ "-^010 ^ 

r<iw _ ^101 _ ^011 _ 1 

•-^110 — ^101 — "-^Oll ~ ^' 

others — 

For double semion model, we have do = 1, di 



(B2) 
-1 and 



^000 _ 1 
"-^000 ^ 

r<ooo _ r< 
'^111 ^ '^i 



001 



/^llO _ ^101 
"-^110 — ^101 

others = 



_ ,-fOll 

— '^100 
_ ,-(011 

— "-^011 



_ ,-.101 _ 

— "-^oio — 

= -1, 



(B3) 



argument of section HI 



Our convention for ordering the indices of G is that the 3 
upper indices i, j, m live on one of the faces of the tetra- 
hedron, while the corresponding lower indices k,l,n live 
on the opposite links (Fig. 12 1. There is no further ambi- 
guity in the index ordering since the G-symbols for these 
models have full tetrahedral symmetry. 

An important property of these partition functions is 
that they are independent of the choice of triangula- 
tion, and depend only on the topology of the space-time 
manifold M. This triangulation independence should 
not be taken for granted: it only comes about because 
^fcin > ^ satisfy highly nontrivial algebraic relations.!^ 

We are now ready to discuss the space-time descrip- 
tion of the duality between the spin models Hq,Hi and 
the string models -ffj.c, H^.s- By analogy with the Hamil- 
tonian description (section |IV| , we place the dual Ising 
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FIG. 12. (Color online) In the Turaev-Viro model ( |Bl| ), 
the degrees of freedom i,j,k = 0, 1 live on the links of the 
triangulation and G^^^ gives a weight to each tetrahedron 
in the space-time triangulation. In the dual spin model, the 
degrees of freedom are Ising spins living on the vertices of the 
triangulation while the dual weight G is defined by mapping 
domain walls between spins onto the link variables i — 0, 1. 
For example, the above configuration corresponds to G(^, t,i 
,t) = Gi\q. The thick red lines denote links with i — I. 



spins on the vertices of the tetrahedra and then define the 
duality by mapping each Ising spin configuration onto its 
corresponding domain wall configuration (Fig. 12). To 



be precise, given any Ising spin configuration, we define 
a corresponding configuration of i = 0, 1 by placing i — 1 
on the links where the adjoining spins are anti-parallel 
(i.e. where there is a domain wall) and z = on the 
links where the spins are parallel (i.e. where there is 
no domain wall). Importantly, we can see from ( |B2|B3P 
that G vanishes for configurations of i,j,k... which do 
not correspond to valid domain wall configurations, so 
this correspondence is (locally) one-to-one. In this way, 
we can map the two Turaev-Viro models (Bl) onto two 
spin partition functions 



Z 



1 

2^ 



'[[d{go,gi) Yi G{go, 91,92, 93), 

ip253--- link tetrahedron 

(B4) 

where 91 =t,i runs over the two Ising spin states. The 
dual G-symbols G{9q, 91, 92, 93) are defined by 
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while d{go,gi) is given by 



d{t,t)=d{i,i)=do 

d(t,i) = d(;,t) = rfi' 



(B5) 



(B6) 



By construction, G and d are invariant under Ising 
symmetry: 

0(990,991,992,993) = 0(90,91,92,93), 

d(g9o ,991) = d(9o,gi) (B7) 

where g- 't=i, 9' i—t- Thus, the dual partition functions 
(B4) both describe Ising symmetric phases. We expect 



that these two phases correspond to the two types of 
paramagnets, Hq,Hi. 

In addition, we note that the dual partition functions 



( B4 ) satisfy the property that they are independent of the 
choice of triangulation: this result follows from the corre- 
sponding property of the Turaev-Viro partition function. 
This property suggests that the two actions described by 
( B4 ) can be regarded as Z2 topological non- linear sigma 



models similar to those constructed in Ref. [11] 

We expect that this construction of topological non- 
linear sigma models can be generalized to arbitrary 
bosonic 2D SPT phases with finite unitary symmetry 
group G. The first step is to find all the Turaev-Viro 
models associated with the group G. In these Turaev- 
Viro models - which are equivalent to the topological 
gauge theories discussed in Ref. |23] - the link labels i 
run over the different group elements of G. (One can 
also construct Turaev-Viro models by placing irreducible 
representations of G on the links, but this approach is 
less convenient here) . One can then construct dual mod- 
els by placing group elements ( "spins" ) on the vertices of 
the tetrahedra, and mapping the domain walls between 
these generalized spins onto the group elements living on 
the links. The result will be a set of topological non-linear 
sigma models with symmetry group G. 



Appendix C: Parity of Ua,Ub under 5* 

In this section, we show that the operators Ua, Ub de- 
fined by UaUbl'^') = can always be chosen so that 
Ua, Ub are either both even or both odd under S. Fur- 
thermore, we show that this even or odd parity must be 
the same for all pairs of endpoints a, b. 

We begin with the first claim - showing that Ua,Ub 
can always be chosen so that they are either both even or 
both odd under S. To derive this fact, we define even and 
odd combinations Ua± — ^(Ua± S^^UaS), and similarly 
for Ub- We then have 

I*) = (Ua+ + Ua-)(Ub+ + Ub-)m (CI) 

We note that = W/j]^) have the same parity 

under S since is even under S. It then follows from 
symmetry that 



Ua+Ub-W) +Ua-Ub+W) 



(C2) 



On the other hand, it is easy to see that the two states 
Ua+Ub-\^') ,Ua-Ub+\^') must be orthogonal to one an- 
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other: 

■ (vl/'lC/^C/b+l*') 
= (C3) 

where the second equahty follows from the fact that \^') 
has short range correlations, and the last equality follows 
from the fact that \'i>') has a definite parity under S. 
Given that the two states are orthogonal and sum to 
zero, both states must vanish: 

Ua+Ub^m = t/a-C/b+l*') = (C4) 

Next, we use the fact that |^') has short range correla- 
tions to deduce that 

= {^'\{Ul^Ua+)iUlU,^)\^') 



Hence, either Ua+\^'} = or ?76_|*') = 0. In the 
same way, we can show that either Ua-\'^') = or 
Ub+l'i'') = 0. There are only two consistent possibilities: 
either Ua^W) = f/b-l*') = or U^+m = C/fc+l*') = 0. 
In the first case, we can replace Ua — > Ua+, Ub Ub+ so 
that both operators are even under S. Similarly, in the 
second case, we can replace Ua — > Ua-,Ub — Ub- so 
that both operators are odd under S. This establishes 
the first claim. 

We now prove the second claim - i.e. that this even or 
odd parity is the same for all endpoints o, b. To see this, 
let /3 be a path joining a, b and /?' be a path joining a with 
some other point c. Then, we have UaUbWpl"^) — \'^) 
and U'JJ'^Wp' \^) = \^) for some operators Ua,Ub,U^,U^. 
Now, by the result above, we know that Ua, Ub have the 
same parity and Ua , U^ have the same parity. Also, it is 
not hard to see that Ua, Ua have the same parity: in fact, 
we can always choose Ua = Ua up to a phase factor, since 
l^*), ly^l^E') have the same local expectation values 
near a. It then follows that Ub, U^ also have the same 
parity under S. 

We now repeat this argument for a path /?" connect- 
ing c with some other point d, letting U'jU'jWfii'l'i!) = 
1^'). By the same reasoning, we find that Ua, Ub, /7", U'J 
all have the same parity under S for arbitrary points 
a, b, c, d. This establishes the claim. 
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